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ABSTRACT 

We investigate microlensing amplification of faint stars in the dense stellar cluster 
in the Galactic Center by the super-massive black hole, which is thought to coincide 
with the radio source SgrA* . Such amplification events would appear very close to the 
position of SgrA* , and could be observed, in principle, during the monitoring of stellar 
proper motions in the Galactic Center. 

We use observations of the near-infrared iT-band (2.2 ^m) luminosity function in 
the Galactic Center and in Baade's Window, as well as stellar population synthesis 
computations, to construct empirical and theoretical K luminosity function models for 
the inner 300 pc of the Galaxy. These, together with the observed dynamical properties 
of the central cluster and inner bulge, are used to compute the rates of microlensing 
events which amplify stars with different intrinsic luminosities above specified detection 
thresholds. 

We present computations of the lensing rates as functions of the event durations, 
which range from several weeks to a few years, for detection thresholds ranging from 
Kq = 16 m to 19 m . We find that events with durations shorter than a few months 
dominate the lensing rate because of the very high stellar densities and velocities very 
near the black hole, where the effective lens size is small. For the current detection 
limit of Kq = 17 m , the total microlensing rate is 3 x 10~ 3 yr _1 . The rate of events with 
durations > 1 yr is 1 x 10~ 3 yr" 1 . The median value of the peak amplification for short 
events is AK ~ 0.75 m above the detection threshold, and is only weakly dependent on 
Kq. Long events are rarer, and are associated with more distant stars, stars at the low 
velocity tail of the velocity distribution, or stars that cross closer to the line-of-sight to 
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SgrA* . Therefore, the median peak amplifications of long events are larger and attain 
values AK ~ 1.5 m above the threshold. 

Recent proper motion studies of stars in the Galactic Center have revealed the 
possible presence of one or two variable A-band sources very close to, or coincident 
with, the position of SgrA* (Eckart & Genzel 1996; Genzel et al. 1997; Ghez et al. 
1998). These sources may have attained peak brightnesses of K ~ 15 m , about 1.5-2 m 
above the observational detection limits, and appear to have varied on a timescale of 
~ 1 yr. This behavior is consistent with long-duration microlensing amplification of 
faint stars by the central black hole. However, we estimate that the probability that a 
single such an event could have been detected during the course of the recent proper 
motion monitoring campaigns is ~ 0.5%. A ten-fold improvement in the detection 
limit and 10 yr of monthly monitoring could increase the total detection probability to 
~ 20%. 



Subject headings: Galaxy: center — Galaxy: kinematics and dynamics — Galaxy: 
stellar content — gravitational lensing — infrared: stars 



1. Introduction 



Recent proper motion studies of infrared-luminous stars in the Galactic Center (GC) ( [Eckart 



fc Genzel 1996] |Genzel et al. 1997| , Ghez et al. 19981 ) have convincingly demonstrated the existence 
of a compact ~ 2.6 1O 6 M0 dark mass in the dynamical center of the GC, which is located within 
0.1" of the radio source SgrA* (|Ghez et al. 1998|) (1" = 0.039 pc at 8kpc). Lower bounds on the 
compactness of this mass concentration, together with dynamical considerations, argue against 
the possibility of a massive cluster, and point towards a super-massive black hole as the likely 
alternative ( penzel et al. 19*97 ; Maoz 1998). This conclusion has important implications for basic 
issues such as the prevalence of massive black holes in the nuclei of normal galaxies, and the 
nature of the accretion mechanism that makes SgrA* so much fainter than typical active galactic 
nuclei QMelia 1994 |Narayan, Yi fe Mahdevan 1995Q . 

Wardle & Yusef-Zadeh (1992) considered several gravitational lensing effects that could 
be induced by a massive black hole in the GC. In particular, Wardle & Yusef-Zadeh pointed 
out that gravitational lensing would occasionally lead to the amplification and splitting of the 
stellar images of stars which happen to move behind the black hole, and that such events would 
typically last from months to years, depending on the distance of such stars from the black hole. 
Wardle & Yusef-Zadeh also estimated that an optical depth of order unity for such events requires 
an observed central surface density > 1000 stars arcsec -2 , and that this in turn would require 
an angular resolution of ^ 0.01" to individually separate the lensed images from the crowded 
background of faint stars. However, the photometric sensitivities and spatial resolutions required 
for such observations are far beyond current observational capabilities. Presently, the deepest 
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near-infrared K-band (2.2 fim) images of the central arcsec reach down to K = 17 m flGhez et 



al. 1998). As we will show, at this limiting magnitude the expected central surface density is 



20 stars arcsec 2 ( Davidge et al. 1997| ). The highest spatial resolutions obtained so far are the 



diffraction limited resolutions of ~ 0.15" ( |Eckart et al. 1995| ; pavidge et al. 199~7| ) and ~ 0.05" 
(Ghcz et al. 1998) achieved in the proper motion surveys. 



In this paper we investigate a different possibility, namely microlensing amplification of faint 
stars by the central black hole. Such events occur when the amplified but unresolved images of 
faint stars rise above the detection threshold and then fade again as such stars move behind the 
black hole close to the line-of-sight to SgrA* . We present a quantitative study of such microlensing 
events, and we compute in detail the microlensing rates as functions of the event durations, 
for a wide range of assumed detection thresholds. In addition, we also consider the possible 
amplification of sources which lie above the detection thresholds, and we re-examine the limit in 
which the two lensed images can actually be resolved. 

Our study is motivated by several recent developments. Deep infrared star counts in the 



inner GC (Blum, Sellgren & DePoy 1996; Davidge et al. 1997) and infrared and optical star 



counts in Baade's Window ( [Tiede et al. 1995 ; Holtzman et al. 1998) now make it possible to 



reliably model the infrared stellar luminosity function in the vicinity of the Galactic Center. The 
on-going proper motion monitoring campaigns of the inner few arcsec of the GC record both the 
positions and fluxes of individual stars in the field at a sampling rate of 1-2 observing runs per 
year. As a by-product, these measurements can be used to search for microlensing events, albeit 
at a low sampling rate. Such events would appear as time varying sources very close to SgrA* . 
It is therefore intriguing that one or two variable IR sources may have already been detected 



close to, or coincident with, the position of SgrA* (Genzel et al. 1997; Ghez et al. 1998). These 
sources may have brightened to K ~ 15 m before fading from view, and appear to have varied on 
a timescale of ~ 1 yr. As we will argue, this behavior is consistent with the expected behavior 
of bright long-duration microlensing events. However, we will also argue that the probability 
that a single such event could have been detected during the course of the recent proper-motion 
monitoring campaigns is small (~ 0.5%). The detection probabilities will increase considerably as 
the observational sensitivities improve (e.g. with the advent of adaptive optics), and if SgrA* is 
monitored more frequently than has been done so far. 

The structure of our paper is as follows. In §^ we set up the formalism required for our 
computations. Some of the more technical aspects are discussed in the appendix. In §|| we discuss 
our treatment of the stellar densities, velocity fields, and ET-band luminosity functions, which 
enter into the computation of the lensing rates. We present our results in §||, where we provide 
computations of the lensing rates as functions of the event durations for a wide range of assumed 
detection thresholds. In §|5| we present a discussion and summary. 
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2. Lensing by the super-massive black hole 

The effective size of a gravitational lens at the lens plane is set by the Einstein radius, R E , 

R, = f 4G f* ^^ V 72 ~2.210 15 (M 2 , 6 d 1 )^cm, (1) 



V c 2 D + d 

where G is the gravitational constant, c is the speed of light, M. is the lens mass (here the black 
hole mass), D and d are the observer-lens and lens-source distances, respectively (see e.g. review 
by Bartelmann & Narayan 1998). We assume, as will be justified below, that d <C D, and define 
M2.6 = M./2.61O 6 M and d\ = d/lpc. The effective size of the lens at the source plane is 
R s = R E (D + d) I D ~ R E . The angular size of the Einstein radius, 8 E , is 

9 E ~ 0.018" Dg^AW!) 1 / 2 , (2) 

where 8-Dskpc is the Sun's galactocentric distance ( parney et al. 1995| ). 

In our study we assume that any IR-extinction associated with an accretion disk or an 
'atmosphere' near the event horizon is negligible. 

We begin by showing that in our problem the stars can be treated as point sources and that 
the linear approximation (small light-bending angle approximation) holds. First, when d <C D, a 
star can be approximated as a point source as long as its radius, R±, is much smaller than R E /A, 
where A is the required amplification to observe the lensed source above the detection threshold 
(see §[2,1|). As we show below (§ |3.2.1 ), the faintest stars that contribute significantly to the lensing 



are low-mass stars with R+ <^ 10 11 cm, which require amplification factors of A <^ 100. Brighter 
stars require progressively smaller amplification factors, down to A ~ 1, which for detection 
thresholds K ~ 17 m correspond to stars with R+ < 10 12 cm. Assuming a mean stellar mass of 
~ 1M©, a central mass density of p ~ 4 10 6 M pc~ 3 in the GC flGenzel et al. 1996Q implies 



a mean stellar separation of 5 ~ 0.005 pc. Even at such a small distance from the black hole, 
R E ~ 1.5 10 14 cm, so that the condition R+ <C R E /A is generally satisfied for all relevant stars in 
the system. We thus conclude that the point source approximation holds over the entire range of 
interest. 

Second, the linear approximation holds as long as the Einstein radius is much larger than the 
Schwarzschild radius of the black hole, R„ 

R E /R. ~c(d/GM) 1/2 ~ 2.810 3 (di/M2. 6 ) 1/2 , (3) 

which even for d = S is as high as ~ 200. 

A point lens produces two images, one within and one outside the Einstein radius. The 
angular separation between the two images is 



A9 = 6 E Vu 2 + 4, 



(4) 
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and their mean angular offset from the optical axis is u9 E , where u9 E is the angular separation 
between the unlensed source and the lens. As will be shown in §||, u -C 1 for amplification above 
present-day detection thresholds in the GC, so that AO ~ 29 E . 

Three angular scales in the problem determine the way the lensing will appear to the observer: 
the Einstein angle 9 E (d), the FWHM angular resolution of the observations, <fi, and the mean 
projected angular separation between the observed stars, Ap(Ko), where Kq is the detection 
threshold. 

The lensed images have to be detected against the background of a very dense stellar system. 
This background will be low as long as as the central surface density of observed stars, So = Ap~ 2 , 
is small enough so that tt^Sq < 1. Thus, for a given detection threshold, an angular resolution of 
at least 

( j ) = Ap(K )/^. (5) 
is required to detect the lensed star. Lower resolutions correspond to the regime of 'pixel lensing' 



( Crotts 1992 ), which we do not consider here. 

When 6 E < (f>, the two images will not be resolved^, and the lensed source will appear as a 
microlensing event. For a given angular resolution, there is a maximal lens-source distance for 
microlensing, 

d »=^tX^) pc - (6) 

More distant stars will have 9 E > 4> and their two images will be separately resolved. 



2.1. The microlensing rate 

The unresolved images of a faint star at d < close to the line of sight will appear as a 
'new source' at the position of SgrA* , to within the observational resolution. The total source 



amplification, A, is related to u by ( Paczynski 1986 ) 



.2 



U 



2A/VA 2 -1-2, (7) 



and for small u, A ~ 1/u. The maximal amplification along the star's trajectory is reached when 
its projected position is closest to the lens. A star of stellar type s and an absolute K magnitude 
K s , which is observed on a line of sight with an extinction coefficient Ak and distance modulus A, 
will be detected above the flux threshold K if it is amplified by at least A s = io-0-4(Ka-K s -A-A K ) _ 
This corresponds to a maximal impact parameter of uq s . As the minimal required amplification 



4 We are assuming here that, generally, a separation of at least 2(f> is required to resolve the two microlensed 
images. The exact value of the minimal separation probably lies in the range to 20, and depends on the details of 
the procedure for faint source separation in the crowded inner l" of the GC. 
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approaches 1, uq diverges. This, and the fact that in practice the threshold is not sharply defined, 
leads us to introduce a cutoff at uq = 1, which implies a minimal amplification of A = 1.34. 

The basic quantity that is required for predicting the lensing properties is the differential 
lensing rate as function of stellar type, event duration, amplification, and source distance from the 
black hole. We derive an explicit expression for the differential lensing rate in the appendix. Here 
it is instructive to discuss the relations between these properties by considering the simpler case of 
the total lensing rate, regardless of the event duration. The total integrated rate for microlensing 
amplification of background stars above the detection threshold is 

r(l?o) = 2(uo> (* R&vin+dT, (8) 
J n 

where r\ is the inner radius of the central stellar cluster, r2 the maximal radius for producing 
unresolved lensed images, V2 is the transverse 2D stellar velocity averaged over the velocity 
distribution function, and n* is the total number density of stars. Here and below, the notation 
(. . .) designates the average of the bracketed quantity over the stellar types with K s > Kq, 
weighted by their fraction in the stellar population, f s , where it is assumed that f s does not 
depend on r. The properties of the stellar population enter the integrated rate only through the 
mean impact parameter (no), which for uq <C 1 is simply 

{s\K a >K } 

~ J2 f s A- l = (F K )/F Q , (9) 

{s\K s >Ko} 

where Fr is the observed (dust extinguished) stellar iT-band flux and Fq is the detection threshold 
flux corresponding to Kq. 

We characterize the microlensing time-scale for stars of type s as the average time they spend 
above the detection threshold, 

f(K ) = ^uq iS R s v 2 , (10) 
where a 7r/4 factor comes from averaging over all impact parameters with u < uq )S . 

The lensing rate, amplification and event duration are inter-related. For uq -C 1, the 
median value of the distribution of maximal amplifications is simply A{uq/2) ~ 2A(uq). Note, 
however, that the mean maximal amplification, (A max ) ~ ^"(uqu) -1 du, diverges logarithmically. 
Generally, a fraction x of the events will have a maximal amplification of Aq/x or more (i.e. a 
peak magnitude above the threshold of AK = 2.5 log x or less). The rate of such events is smaller, 

r(x) = xT(Kq) . (ii) 

The time-scale of events amplified by a factor of 1 jx above the threshold is somewhat longer than 
the average time-scale (Eq. 10), since they must cross closer to the line of sight, 

2 1 I k sin _1 x\ , „ . , , 

f{x) = - Vl-x 2 + f(K ) , (12) 
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which approaches At(Kq)/it for large amplifications. Conversely, for a given V2 and uo, s , even a 
small increase in f is associated with a considerable increase in the maximal amplification. 



2.2. Resolved lensed images 

When 9 E > (f> at d > d^, the two images can be resolved. As we show in §||, d^ is already 
quite small for present-day angular resolutions, and will become smaller still as the resolution 
improves. This implies that there may be a non- negligible contribution to the lensing rate from 
regions beyond cL. We therefore have to consider also the case of resolved images. 

For aCl, the two images will appear at an offset of ~ 9 E from SgrA* , on opposite sides of it. 
The amplifications of the individual images are related by A = A\ + A2 and A-i = A\ + 1, which 
for the high amplifications that are requires in the GC can be approximated as A\ ~ A2 = A/2. 
The formalism used for calculating the lensing of unresolved images can therefore be applied 
in this case simply by raising the effective detection threshold by a factor of two (0.75 m ). This 
makes resolved images harder to detect, but on the other hand, if both images are observed, the 
identification of the event as lensing is much more certain. 



2.3. Lensing of observed bright stars 

The formalism for describing unresolved and resolved lensing of stars from below the detection 
threshold can be also extended to cases where the microlensed source is an observed bright, 
K < Kq star. At the current detection threshold, the observed central surface density is ~ 10 
arcsec -2 flGenzel et al. 1997 ; Ghez et al. 1998| ). For such a small number of stars, whose orbits 



can be tracked individually, a statistical treatment is not very useful. However, our calculations 
indicate that a ten-fold improvement in the detection threshold will yield an observed central 
surface density of at least ~ 100 arcsec -2 (Eq. |25| ). For such a large sample, a statistical treatment 
is more meaningful. We therefore present below, for completeness, results for microlensing 
amplification of stars fainter than the detection threshold ('faint-star lensing') as well as for stars 
brighter than the detection threshold ('bright-star lensing'). For bright-star microlensing we 
assume that the effective minimum magnification factor for detection is A = 1.34 (uq = 1). 



3. Modeling the stellar population in the Galactic Center 

Three basic quantities enter into the computation of the microlensing rate: the stellar number 
density distribution, the stellar velocity field, and the if -band luminosity function (KLF). 

The stellar population in the central ~ 100 pc appears to consist of a mixture of old bulge 
stars, and 'central cluster' stars which may have been produced in various star-formation episodes 
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during the lifetime of the Galaxy (Genzel, Hollenbach & Townes 1994; |Serabyn fc Morris 1996| ). 



Evidence for recent star formation in the GC has been provided by Krabbe et al. (1995) who 
found a concentration of luminous early-type stars within a few arcsec of SgrA* , implying a recent 

10 7 yr) starburst in which ~ 10 3,5 stars were formed. Additional young stellar systems, the 
'Arches' and 'Quintuplet' clusters also exist close the GC and contain large numbers of massive 
stars ( [Serabyn, Shupe fc Figer 1998] ). A further indication of on-going star- formation in the GC 
is the fact that the KLF in the central cluster extends to more luminous stars than the KLF of 
the Galactic bulge as observed via Baade's window ( [Lebofsky Sz Rieke 1987| ; [Tiede et al. 1995 ; 



Blum, Sellgren &: DePoy 199€ ; Davidge et al. 1997|) . Serabyn & Morris (1996) suggested that 



continuous star formation in the central cluster is maintained by molecular clouds in the GC, and 
that the ~ r~ 2 radial light profile of the central cluster reflects the distribution of the star-forming 
molecular clouds. 

In our analysis we make the simplifying assumption that the discrete young clusters in 
the GC can be modeled by a volume averaged and smoothly distributed stellar population. In 
particular, we do not consider here lensing events that might be associated with the innermost 
stellar cusp, which has been identified by Eckart & Genzel (1997) and Ghez et al. (1998) in the 
immediate vicinity of SgrA* . An over-density of stars above the smoothed distribution very near 
the black hole may contribute to very short duration microlensing events. However, the total 
stellar mass and luminosity function of the stars that are associated with this 'SgrA* cusp' are 
poorly constrained at the present time, so that computations of the lensing rates require more 
extensive modeling and analysis. We will present such computations elsewhere. 

Because the lensing time-scales increase with distance, long duration events, which are 
relevant for low sampling rates of the current proper motion studies, tend to be associated with 
more distant stars. We therefore consider both stars in the central star-forming cluster, as well 
as more distant old-population bulge stars in our analysis. We now proceed to discuss the stellar 
densities, velocities and KLFs of these two components. 



3.1. Stellar densities and velocities 

Genzel et al. (1996, 1997) derived density and velocity models for the central cluster based 
on fits to the observed star counts, stellar radial velocities and proper motions in the inner few pc. 
Their mass density distribution is parameterized by a softened isothermal distribution 

Pcore(r) = - P ,° . . 2 (13) 
1 + 3(r/r c ) / 

Where r c is the core radius and p c is the central density, with best fit values of r c = 0.38 pc 
and p c = 410 6 M Q pc~ 3 . The ID velocity dispersion of the late type stars in the core, which are 
dynamically relaxed, is modeled by ( penzel et al. 19"96| ) 

<W = (55 2 + 103 2 (r/r 10 )- 2o ) 1/2 kms" 1 , (14) 
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where a = 0.6 and no is the projected distance corresponding to 10". 

In the inner GC, the mass density is strongly constrained by the dynamics, whereas the K 
luminosity density v (Lkq pc -3 ) is not well defined due to the patchy nature of the extinction. 
The situation is reversed on scales larger than 1 kpc, where the observed rotation curve and 
velocity dispersion are harder to interpret, but u can be de-projected from the observed surface 
brightness. Kent (1992) proposed that both the rotation curve and the surface brightness along 
the major axis of the Galactic bulge in the mid-plane of the galaxy at r ^ lkpc can be described 
as a superposition of bulge and disk components with a mass-to-light ratio[] T = 1T Q and a 
luminosity density 

v{r) = i-Wge + ^disk (15) 
= 3.53K (Jr) + 3exp (-^j L KQ p C - 3 , 

where here Kq is a modified Bessel function (not to be confused with the detection threshold), 
r& = 667 pc and = 3001 pc. Kent also suggested a v oc r -1 ' 85 extrapolation of v towards 
the center based on observations of the intensity variation in the inner 10 pc. We replace this 
extrapolation with an updated and non-divergent model by adopting the Genzel et al. (1996) 
mass density model (Eq. |i~3| ) for the core and by assuming T = 1T for the bulge and disk, so that 
Pbuige = i^buige and pdisk = ^disk- We further assume that the bulge is axisymmetric, and model the 
mass density over the entire distance range by 

P = Pcorc + Pbulge + Pdisk • (16) 

Figure [i] shows our mass density model, and in particular the emergence of the bulge component 
at ~ 100 pc. The central cluster completely dominates the total mass density in the inner 10 pc, 
but then falls to ~ 85% of the total at 50 pc and to ~ 65% at 100 pc. 

The velocity field on this scale includes both the bulk rotation and the random motion. 
We approximate Kent's models for the galactic rotation and the velocity dispersion in the inner 
300 pc along the major axis of the bulge by log-linear fits. We assume that the rotation velocity is 
perpendicular to the line of sight, and that its contribution to the transverse velocity is 

v mt = 80 + 201og 10 (r/lpc) kms" 1 , (17) 

and that the ID velocity dispersion in the bulge is 

o- bulge = 60.9 + 18.91og 10 (r/lpc) krns" 1 . (18) 



5 Following Kent (1992), we define the mass-to-light ratio as T = (M /Mq) / (Lk / Lkq), where the solar 
monochromatic luminosity at 2.2/im is Lkq = 2.154 10 32 ergs -1 /im" 1 . Note that Genzel et al. (1996) define 
this quantity as T' = (M /Mq)/(XL\/ Lq), with A = 2.2/im. The two definitions are related by T' = 8.07T. 
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We model the ID velocity dispersion over the entire distance range by 

a = max(o- cor c, cr bulge ) • (19) 

Note that both the proper motion of SgrA* , which is <J 20kms~ 1 (Backer 1996), and the 



Sun's galactocentric rotation have a negligible contribution to the relative proper motion of the 
source and lens because of the high stellar velocities near the dynamic center and because d <C D. 



3.2. K luminosity function 

In our analysis we consider a model for the KLF which is based on the observed portions 
of the KLFs in the central cluster and in the bulge, and on a theoretical KLF computed in a 
population synthesis model for the central cluster. 



3.2.1. The power-law KLF 

The KLF of the bulge has been observed through Baade's window, which at I = 1.0°, 
b = —3.9°, is ~ 0.6 kpc from the GC at the tangential point. The bulge KLF approximately follows 
a single power-law dlog N+fdlog Lk = /3, with (3 = 1.695 ( Tiede et al. 199~5| ), from Mk ~ — 7.5 r 



-m 



to Mk ~ 2 m [], assuming a distance modulus A = 14.5 and Ak ~ 0. The faintest observed stars 
are at the detection threshold, and the power-law KLF likely extends, or steepens slightly, to lower 
luminosity stars near the main-sequence turn-off point. 



Recent HST observations of the optical luminosity function in Baade's Window ( Holtzmarj 



et al. 1998) probe it down to very low stellar masses (~ 0.3 M ), well below the main-sequence 



turn-off point. The V-band luminosity function (VLF) presented by Holtzman et al. (1998) shows 
that a sharp break occurs at the turn-off point My ~ 4.5-5 m , corresponding to ~ 1 Mq stars. We 
extrapolated the power-law KLF of Tiede et al. (1995) down to Mk = 3.5 m (assuming V-K ~ 1.5 
for IMq stars) and compared it to the Holtzman et al. (1998) VLF at My = 5 m . We find that 
the predicted K star-counts somewhat underestimate the V star-counts, but agree with them to 
within a factor of 2. We thus conclude that the power-law KLF can be extended down to at least 
Mk = 3.5 m . The observed VLF can be used to determine the behavior of the KLF at even lower 
luminosities by using the V-to-K conversions for low mass stars presented by Henry & McCarthy 
(1993). Doing this shows that the KLF likely flattens at a break point close to Mk = 3.5 m , with 
(3 ~ 1.5 in the range 3.5 m <^ Mk <^ 6.5 m , and turns over at Mk ~ 7 m . 

As we have discussed above, the stellar population in the central cluster appears to be 
consistent with continuous star-formation throughout the Galaxy lifetime (see e.g. Serabyn &; 



s Mk = -2.5 logL if + 84.245 mag for Lk in ergs 1 /im 1 . 
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Morris 1996). However, despite the many differences between the populations in the bulge and 
the central cluster, the observed KLF in the central 12.4" x 11.9" of the GC follows a power-law 
similar to that of the bulge KLF ( pavidge et al. 1997 ), but extends to more luminous stars ( [Blum 



Sellgren fc DePoy 1996 ). The exact upper luminosity cutoff is not well-defined due to statistical 
fluctuations in the counts and the contribution of asymptotic giant branch (AGB) stars, which do 
not follow the power-law. 

The central cluster's KLF is not known beyond the current detection threshold of Kq ~ 17 m , 
which is the region of interest for the microlensing calculations. However, the similarity of the 
KLF power-law indices in the central cluster and bulge KLFs suggests that the two KLFs are 
similar, apart for having different upper luminosity cutoffs. The much lower extinction in Baade's 
window, which is only Ax = 0.14 m as compared to 3.4 m in the GC ( [Rieke, Rieke fc Paul 1989| ), 
makes it possible to extend the KLF in the GC to lower luminosities. Indeed, Blum et al. (1996) 
find that the bulge KLF, after correcting for the Ak difference, can be smoothly joined to the GC 
LF in the central 2' x 2'. Their resulting composite LF has a best fit power-law index j3 = 1.875, 
and extends from Mr ~ — 10 m down to Mk = 2 m . Our theoretical model for the central cluster, 
which we discuss below, indicates that this power-law character is further maintained down to 
M K ~ 3.5 m . 

We therefore adopt the same (3 = 1.875 power-law index for both the central cluster and bulge 
KLFs. For the central cluster KLF we set the upper luminosity cutoff at Kh = 8 m , and for the 
bulge KLF at = 10. 5 m . We set an effective low luminosity cutoff for both at K\ = 21. 5 m , which 
corresponds to stars with masses ~ 1M@. As we show in § |3,2.3 , the lensing rates are insensitive 
to the exact choice of the low luminosity cutoff. 

We note that the observed bulge KLF is better fitted by a somewhat flatter power-law than 
(3 = 1.875, which implies less faint-star lensing candidates. On the other hand, the comparison 
with the optical luminosity function suggests that the KLF is a conservative estimate of the 
total number of stars. In addition, the flatter power-law fails to account for the strong excess 
of horizontal branch (HB) stars above the power- law ( [Tiede et al. 1995| ). These are important 



potential microlensing sources, as they lie in a magnitude range that is just below the threshold if 
the bulge population is observed through the GC extinction of Ak = 3.4 m . We therefore consider 
these small discrepancies in slope and number to be within the limitations of the power-law 
approximation and the observational uncertainties. 



3.2.2. The theoretical KLF 

As a check on the empirically based pure power-law KLF, we have also computed a series of 
theoretical KLFs using our 'population synthesis' code (see e.g. Sternberg 1998|) . 



In our models we used the Geneva stellar evolutionary tracks ( [Bchaerer et al. 1993| ), and 
concentrated on stellar models with twice-solar metallicities, as is indicated by the enhanced 
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abundances of the gas in the GC ( [Morris Serabyn 1996] ). However, recent measurements 



of stellar spectra of cool luminous stars in the GC point to solar metallicities (Ramirez et al 



1998 ). We therefore considered also such stellar models, and verified that our results do not 
depend strongly on the assumed metallicity. We computed the stellar K luminosities using the 
empirical bolometric corrections and V-K colors for dwarfs, giants and supergiants compiled by 
Schmidt-Kaler (1982), and Tokunaga (1998). The Geneva tracks for intermediate mass stars 
(~ 2-7 'Mq) do not extend beyond the end of the early AGB phase. We extended these tracks 
to include also the more luminous thermal-pulsing AGB phases following prescriptions described 
by Chariot & Bruzual (1991) and Bedijn (1988). In our models we assume explicitly that stellar 
remnants and stars less massive than 0.8M & are negligible sources of iT-band luminosity. This 
low-mass luminosity cutoff roughly corresponds to the low luminosity cutoff of the empirically 
based power-law KLF. 

We constructed theoretical models for a range of cluster parameters assuming continuous 
star-formation lasting for 10 Gyr. We considered models with pure power-law IMFs, and the 
Miller-Scalo IMF ( [Miller fc Scalo 1979 ; Scalo 1986| ) for a range of lower- and upper-mass cutoffs. 



We selected the model which yields a KLF which best matches the observed KLF for the central 
cluster measured by Davidge et al. (1997). 

We find that a model with a Miller-Scalo IMF ranging from 0.1 to 120 Mq provides the best 
overall fit to the data. In Figure [2] we compare our model KLF for the central cluster with the 
(3 = 1.875 power-law fit of Blum et al. (1996) to their composite KLF. The overall agreement 
between the theoretical KLF and the power-law KLF is remarkable, for both stellar models 
with solar and twice-solar metallicities. The model successfully reproduces both the power-law 
character and slope of the observed KLF of the central cluster. Furthermore the model shows that 
the power-law KLF likely extends down to at least K = 21.5 m . In this model, the most probable 
lensed sources are K ~ 21 m stars with M* ~ 1M , ~ 10 11 cm just off the main sequence, 
which require a magnification of order A ~ 50 to be detected above the threshold. 

Our population synthesis model predicts a mass-to-light ratio T = O.24T0, in excellent 
agreement with T ~ 0.25 T measured in the central cluster by Genzel et al. (1996). 



3.2.3. Normalizing the KLF 

As shown by eq. pi, the lensing rate depends on, n+, the total number density of stars which 
are effective sources of K-band luminosity. We refer to such stars as 'K-emitting' stars, and define 

n* = 4p, (20) 
m 

where p is the total dynamical mass density, fh is the mean stellar mass of the K-emitters, and 
is the fraction of the total dynamical mass contained in K-emitting stars, i.e. excluding low-mass 
stars and remnants (i.e. objects which lie below the effective low-luminosity cutoff of the KLF) 
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and gas clouds. 

The observed star counts, dN° hs /dLx, within an area S, can be used to obtain a best fit value 
of f*/fh given the constraint 

div° bs df u r 

-77 — = —j—— pdSdr, (21) 
dLx dLx m J 

where the KLF df /dL^ is normalized to 1. The mass-to-light ratio over the integration volume 



used in Eq. 21 can then be deduced from the fit by 



™ M fh 

T = — = TT - ■ (22) 
Lk ULk 

where Lk is the average K-band luminosity of the K-emitting stars. 

The power-law KLF we adopt in our computations (df/dLx oc L^f) is characterized by 
1 < P < 2, Li -C L u and L\ -C Lq <C L u , where L/ and L u are the effective lower and upper 
luminosity cutoffs to the KLF, and Lq is the K luminosity corresponding to the detection 
threshold. By approximating u ~ 1/A (where A is the required amplification) the mean stellar 
luminosity and mean impact parameter for such KLFs are given by the simple approximate 
relations 



Lk 



(2-/3) L )- p 



(2-P) \L 

It also follows that n*, T and the lensing rate V scale as 



1-0 



1 

p-r 

T ~ (2-P)Lt\ 

r ~ ( Uo )n, ~ -^-L^ . (24) 



2-/3 



Eq. 24 shows that the total lensing rate is insensitive to the low-luminosity cutoff of the KLF. 
This simply reflects the fact that an increase in the number of if -emitting stars as L\ decreases 
is offset by a correspondingly smaller mean lensing impact parameter for the stellar system, and 
vice-versa. This important property allows us to robustly compute the lensing rate even if the 
effective lower luminosity cutoff of the power-law KLF is not well determined^]. 

We note that since the differential contribution of stars with luminosity Lk s to the mean 
impact parameter (uq) scales like {Lk, s / L$)L~ff s ~ L^g, the integrated contribution of stars from 



7 We note that if /3 > 2 then the rate does depend on Li with Y ~ Lf~ /3 . However, as we have discussed, the 
observations of Holtzman et al. (1998) strongly suggest that the KLF flattens, rather than steepens, below our 
assumed value for Li. 
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a l m -wide bin is, for f3 = 1.875, a very weakly decreasing function of the bin's K magnitude. We 
therefore expect that the lensed stars will exhibit a wide range of intrinsic K magnitudes, with a 
weak trend towards lensing by stars close to the detection threshold. 

It is uncertain at which radius the stellar population makes the transition from a population 
that is characteristic of a star-forming cluster to a more bulge-like population. However, as we 
argued above, the observed properties of the KLFs in the GC and the bulge suggest that they 
are very similar for K > 10. 5 m . Since the normalization of the KLF does not depend strongly on 
the upper luminosity cutoff (Eq. [24]), and since the very high luminosity tail of the KLF is not 
relevant for the lensing calculations (the 8 m < K < 10. 5 m stars have a negligible contribution to 
the lensing rate of observed stars), we infer the value of f*/fh using the star-counts observed in 
the core, and we adopt this normalization for the entire inner 300 pc. 

In carrying out this procedure^] we used the observed number counts in the K = 14 m bin 
(stars per mag per arcsec -2 ), averaged over the 12.4" x 11.9" field observed by Davidge et al. 
(1997). The star counts in this luminosity range are likely complete, and this range is also well 
separated from the regions that are affected by AGB and HB stars, which cause deviations from 
the power-law behavior (Tiede et al. 1995). Using Eq. [H], we then find that f*/m = 0.2 Mq 1 . 
This value for f*/rh can be reproduced, for example, by a choice of /* = 0.2 and fh = IMq, which 
is comparable with the values /* = 0.22 and m = 0.84M Q of our population synthesis model for 
the central cluster. For the central cluster power-law KLF, Lk = 22Lkq, so that Eq. ^2] yields 
T ~ 0.22 Tq, again in excellent agreement with T ~ 0.25 T inferred by Genzel et al. (1996) for 
the central cluster. 



4. Results 

We have carried out detailed computations of the lensing event rates using the formalism 
described in §|2| and in the appendix, and using the stellar number and velocity distributions and 
the power-law .fT-band luminosity function, as discussed in §|3[ The integrations were carried out 
from a minimum distance r\ = 0.005 pc, equal to the mean central stellar separation, to a distance 
ri = 300 pc, where the integrated lensing rates approach their asymptotic values. 

The normalized central cluster KLF allows us to estimate the stellar surface density as a 
function of the detection threshold, and from it derive the mean angular separation of the stars, 
Ap, the required angular resolution of the observations <fi, and the maximal distance for observing 
(unresolved) microlensed stars (as given by Eqs. |5| and ^). We find that for the central cluster 



The volume integration in Eq. |2l| was carried out by approximating the rectangular field with a circular field of 
projected radius p having the same area. The integration was done over a cylinder of radius p, centered on the black 
hole and extending along the line of sight 300 pc in each direction, a distance large enough to ensure that the surface 
density reaches its asymptotic value. 
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KLF 

logAp = 2.31-0.175^0, 

log So = -4.62 + 0.35^0, (25) 

for Ap in arcsec and So m arcsec -2 . Thus, for a detection threshold of Kq = 17 m we expect a 
central surface density of 21 stars per arcsec -2 for complete counts^]. 

We present our results in a way which makes it possible to flexibly estimate the actual 
detection probabilities for a broad range of observing strategies. We consider the generic 
monitoring campaign to consist of a series of n very short observing runs which are carried out 
during a total time period T (typically T ~ several years), with a mean interval AT between 
each observing run (typically AT ~ a month to a year). We now wish to distinguish between 
'long' and 'short' events, where we define the event duration as the time the source spends above 
the detection threshold. Long events are those with durations r > AT, and would appear as 
time varying sources that brighten and fade during the course of several observing runs. The 
microlensing origin of long events could be verified, in principle, from the shape of the light 
curve and its achromatic behavior. Short events are those with durations r < AT, and would be 
observed as a single 'flash' provided they occur within a time r of any one of the n observing runs. 
In the limit of small event rates the detection probability may be written as 

P -fshort -flong ^TshortFshort -^Tlong, (26) 

where r snor t and Fi ong are the rates of short and long events, and T sn0 rt is the rate averaged lensing 
duration for short events. In the (ideal) limit of continuous monitoring, ri ong approaches the total 
event rate and P = Tr to tal- 

The results of our computations are displayed in Figs. |] and ||. In Fig. || we plot the 
cumulative rates, ri ong (r > AT), for all lensing events with durations r longer than the timescale 
AT, as functions of AT. We present results for events which produce unresolved and resolved 
images for stars which are either intrinsically below or above the detection thresholds. We present 
computations for detection thresholds ranging from 16 to 19 magnitudes. The total lensing rates, 
can be read off the plot from the asymptotic values of the curves as AT — ► 0. The curves are 
flat for timescales less than ~ 1 week, which is shorter than most of the event durations. As AT 
increases ri ong decreases as a smaller fraction of the lensing events satisfy r > AT. As an example, 
Fig. |] shows that for Kq = 17 m , the total lensing rate is equal to 3 x 10~ 3 yr -1 , and that for 
events with durations greater than 1 yr the rate is equal to 1 x 10 -3 yr -1 . 

In Fig. |] we also plot the rate-averaged lensing timescale, T snort , for events with r < AT, 
which is needed to estimate the detection probability of short events. The values of f s h or t are 



9 Genzel et al (1997) and Ghez et al (1998) reported So ~ 20 and 15 stars arcsec -2 , respectively. However, the star 
counts are incomplete close to the detection thresholds, and the measured So actually indicate a significant density 
enhancement (the 'SgrA* cusp') relative to the immediate surroundings. 
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almost independent of Kq, since the shape of the cumulative rate function is insensitive to Kq. 
We note also that for small timescales Tshort ~ AT/2, as would be expected for a rate which is 
independent of the timescale. The rate of short events is simply r s hort(AT) = riong(O) — ri ong (AT). 
Thus, for Kq = 17 m the rate of events lasting less than 1 yr is 2 x 1CP 3 yr -1 , and the average 
duration of such events is about 3 months. 

In Fig. [| we plot AK\ ong , the median value of the maximal amplifications for long events 
(t > AT), as a function of AT. Long events which amplify sub-threshold stars are associated 
with stars at greater distances from the black hole, stars at the low velocity tail of the velocity 
distributions, or stars that cross closer to the line of sight to SgrA* . Because of the latter 
effect, AK\ ong is greater than 0.75 m , which is the median value for all the events. This effect is 
less marked for resolved lensing, which is characterized by longer timescales, and is very weak 
for lensing of stars which lie above the detection threshold, where AK\ ong approaches the limit 
~ 0.75 m . 

Figure || shows the contributions to the lensing rate from different regions in the GC, for a 
specific example where Kq = 16. 5 m and AT = lyr. Several important features are apparent in the 
results displayed in Fig. || and ||. First, it is evident that the cumulative rate r to tai(< r) approaches 
its asymptotic value at r ^ 10 pc, and that short lensing events due to stars near the black hole 
dominate the total rate. For example, Fig. || shows that the median timescale of unresolved 
amplifications of sub-threshold stars for is ~ 2 months, and Fig. || shows that the median distance 
from the black hole is ~ 0.3 pc. The run of the differential rate for long events, dT\ ong /dr, with 
distance illustrates that the inner regions hardly contribute any long events. Second, unresolved 
lensing of sub-threshold stars have the shortest timescales. This is because such events are due 
mainly to stars close to the black hole where the velocities are high and the lensing cross sections 
are small. Unresolved lensing of stars which lie above the threshold, for which the cross sections 
are larger, have somewhat longer timescales, and resolved events, which are due to stars at larger 
distances from the black hole, are longer still. Third, the rates for resolved lensing are about an 
order of magnitude smaller than those for unresolved lensing events. 

We now apply our results to estimate the probability that the variable iT-band source (or 
sources) reported by Genzel et al. (1997) (source S12) and Ghez et al. (1998) (source S3) was 
a microlensing event. The monitoring of proper motions in the GC has been going on for about 
T = 6yr, at a sampling interval of AT = lyr, with a detection threshold of Kq = 16. 5 m and a 
FWHM resolution of 0.15" flGenzel et al. 1997| ). For this detection threshold, Ap = 0.25" (which 
corresponds to a required spatial resolution <fi = 0.15"), and d^ = 65 pc. We can now use Figs. |] 
and |E] and Eq. 26 to estimate the detection probability and typical amplification of a lensing event 



in this experiment. The rate of unresolved and long events which amplify sub-threshold stars is 
7.5 x 10 yr, with a resulting detection probability -Pi on g ~ 0.5%. The median amplifications of 
such events is Ai^i ong ~ 1.5 m above the detection threshold. The probability for detecting a short 
amplification of sub-threshold stars is -P s hort ~ 0.2%. The probabilities for detecting unresolved 
events involving above threshold stars is of the same order of magnitude. The probabilities 
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for detecting long resolved events are an order of magnitude smaller, and the probabilities for 
detecting short resolved events are negligible in this experiment. Thus, we find that the behavior 
of the variable source (or sources) at SgrA* , i.e. a brightening of a previously undetected source to 
1.5 m -2 m above the threshold with an event duration of ~ 1 yr, is the typical behavior that would 
be expected for a microlensing event. However, we also find that the a-posteriori probability for 
detecting such an event is only ~ 0.5%. 

The probability for detecting a lensing event can be increased considerably by carrying 
out more sensitive observations at higher sampling rates. For example, 10 years of monthly 
observations with a detection threshold of Kq = 19 m will require a resolution of <p = 0.06", 
which is already available ( |Ghez et al. 1998| ), and will increase the total detection probability 
of long lensing events to Pi on g ~ 20%. These estimates are uncertain by a factor of few due to 
uncertainties in the stellar density distribution, the K luminosity function and the extinction. 



5. Discussion and summary 

In the past several years, very high resolution, very deep IR observations of the GC have 
regularly monitored stellar motions within few arcsec of the radio source SgrA* . The primary 
objective of these efforts is to dynamically weigh the central dark mass and set lower bounds on 
its compactness. As a by product, these observations can detect and record the light curves of 
faint time-varying sources in the inner GC over timescales of years. These proper motion studies 
have recently revealed the possible presence of one or two variable X-band sources very close to, 



or coincident with, the position of SgrA* ( Genzel et al. 1997 ; Ghez et al. 199 



The first issue to resolve, when considering results from such technologically challenging 
observations, is whether these sources are real, or merely artifacts of the complex procedures for 
obtaining deep diffraction limited images in the crowded GC. Assuming that such a source is 
indeed real, and that it is not simply a variable star that happens to lie very near to the line 
of sight to SgrA* , there are two interesting possibilities to consider, both directly linked to the 
existence of a super-massive compact object in the GC. A time variable IR source, coincident with 
SgrA* , may be the IR counterpart of the radio source, with the IR flare resulting from fluctuations 
in the accretion process, which has up to now eluded detection in any band other than the radio 



(see e.g. Backer 1996). Another possibility is that the new source is a faint star in the dense 
stellar cluster in the GC, which is microlensed by the black hole. Such amplification events would 
appear as time varying sources very close to the position of SgrA* . 

We note that a detection of microlensing could provide an independent probe of the 
compactness of the central dark mass. The innermost observed stars in the GC set an upper 



limit on the size of the dark mass, r^m ~ 710 R,, where R, is the Schwarzschild radius ( Ghez 



et al. 1998| ). Microlensing has the potential, given a well sampled light curve, to improve on 



the dynamical limit, since for typical values of the lens-source distance, <C R E < 3 10 4 R, 
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(see Eq. |3|). If the central mass is not a black hole, but rather an extended object, then the 
microlensing light curve will deviate from that of a point-mass lens when r^m ^ -Re- Because the 
innermost observed stars already constrain 9^ m <J 0.1", only lensed stars in the inner few pc have a 
small enough 9 E to probe possible structure in the dark mass distribution (Eq. ^). Therefore, the 
question whether there is a significant chance to detect such events is related to the yet unresolved 
question of whether there is a strongly peaked stellar cusp in the innermost GC. We will discuss 
these issues elsewhere. 

In this paper we investigated the possibility of microlensing amplification of faint stars in the 
dense stellar cluster in the GC by the super-massive black hole, which is thought to coincide with 
the radio source SgrA* . We calculated in detail the rates, durations and amplifications of such 
events, and considered separately the cases of unresolved and resolved images and of intrinsically 
faint and bright sources. We presented our results in a general way that can be used to estimate 
the detection probabilities of microlensing for a wide range of observing strategies. 

The background stellar surface density increases with the detection threshold Kq. This 
determines the observational angular resolution required to detect a microlensing event, and fixes 
the maximal distance behind the black hole for which the two lensed images of a star will appear 
unresolved. This maximal distance occurs before the integrated lensing rate reaches its asymptotic 
value, even for present-day detection thresholds. We therefore considered two manifestations of the 
lensing: unresolved microlensing of stars near the black hole, and resolved lensing of stars farther 
away. We find that short lensing events, due to stars close to the black hole, dominate the total 
lensing rate. This reflects the fact that the high stellar density and velocities near the black hole 
more than compensate for the smaller lensing cross-sections there. For this reason, and because 
unresolved images are on average twice as bright as the resolved images, unresolved microlensing 
dominates the lensing rate. We have also considered the lensing amplification of bright observed 
stars. The contribution of this type of microlensing to the total rate becomes progressively more 
important as the detection threshold decreases, and at low sampling rates, which are primarily 
sensitive to longer events. We find that low sampling rates significantly bias the detection towards 
high amplitude events. 

Our predicted lensing rates are small, but not so small as to be negligible. In particular, 
longer, deeper proper motion monitoring done at higher rates, e.g. 10 years of monthly monitoring 
with Kq = 19 m , may have a significant chance of detecting such an event. 

Finally, could either of the variable sources reported by Genzel et al. (1997) (source S12) 
and by Ghez et al. (1998) (source S3) be the amplified microlensed image of a faint star? The 
lack of evidence for related variability in the radio and X-rays, as would be expected in some 
accretion scenarios if the K-B.aie were due to fluctuations in the accretion process, argues against 
the possibility that the new source is the IR counterpart of SgrA* . 

Another possibility is that the new sources are variable stars, which are below the detection 
limit in their low state. Long Period Variable stars (LPVs), which are probably luminous Mira 
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variables, are observed in the GC (Haller & Rieke 1989; Tamura et al. 1996). Typical amplitude 



variations of AK ~ 0.15 m -0.5 m are observed over 1-2 yr, although in some cases the variations 
are as large as AK ^ l m . Haller & Rieke (1989) find 12 LPVs in a 4.5 arcmin 2 survey of the GC 
(not including the central V x 1') down to K = 12 m (Mk = — 5.9 m ). Of these, only one exhibited 
high amplitude variations (l m in 4 months). At K ~ 15 m , the new variable K source is much 
fainter than the LPVs observed in this survey. If it is an intrinsically bright LPV, it must be lie on 
a highly extinguished line of sight. Using the observed surface density of Mk < — 5.9 m , AK ~ l m 
LPVs, it is possible to make a rough estimate of the probability for finding such a star within 
0.15" of SgrA* . Even after taking into account a factor ~ 40 difference in the surface mass density 
between the dynamical center and the survey area at ~ 2.5' from the center, the probability is 
only ~ 0.005%. This of course does not rule out the possibility of an intrinsically faint but highly 
variable star. In any case, if the new source is a variable star, future observations should detect 
continued variability from this source. 

We can also argue statistically against the possibility that these flares are due to microlensing 
of a star by another star in the GC, which happens to be close to the line-of-sight to SgrA* . 
Because m/M, ~ 10 -6 , the typical lensing timescale by a star will be of the order of 1 hr, 10 3 times 
shorter than that of lensing by the black hole (Eqs. [I] and ||). Furthermore, it is straightforward 
to show that the total rate of lensing of a star by a star within an angle 0.1" of the line-of-sight to 
SgrA* is ~ 5 10~ 5 yr -1 for Kq = 17 m , which is 100 time smaller than the rate of lensing by the 
black hole. 

Our analysis has shown that the behavior of the variable X-band source (or sources) at SgrA* , 
in particular a brightening of a previously undetected source to ~ 1.5-2 m above the threshold, on 
a timescale of ~ 1 yr, is the typical behavior that would be expected for a microlensing event. 
However, we also estimate that the probability that such an event could have been observed 
during the course of the proper motion studies that have been carried out thus far is only ~ 0.5%. 
While this probability is small, it is not so small as to rule out this possibility entirely. The 
probability of detecting a microlensing event at SgrA* will increase considerably in the future as 
the observational sensitivities and the monitoring sampling rates improve. 
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A. The microlensing rate and amplification of long duration events 

Figure |] shows that the total lensing rate is dominated by short events (shorter than a few 
month). Such events are shorter than the current mean time between observations, and can 
therefore be observed only in the rare cases where they occur simultaneously with an observing 
run (§||). Of more relevance are long events, which span two or more observing runs. Longer 
events occur when the lensed star is either slow, or when it is far away from the black hole, so 
that its trajectory through the Einstein ring is long (Eq. |l|), or when it passes close to the optical 
axis (Eq. |i~2] ). In this appendix we calculate the rate of events that are longer than some specified 
time, and estimate their median maximal amplification. 

Microlensing events with r > AT are those with transverse velocity v 2 and impact parameter 

uR s such that 

2Rjul >s -u 2 

t= >AT, (Al) 

where s is the stellar type and uq jS is the maximal impact parameter for amplification above the 
threshold. The local rate per stellar type is 

" 2R s n* / du *72DF2(v 2 )(fV2, (A2) 



drds 



where f max = 2^Ju 2 )s — u 2 R s /AT and DF2 is the 2D distribution function of velocities. The 2D 
velocity v 2 is composed of an ordered component v rot , which we assume to be perpendicular to 
the line of sight, and a random isotropic component, which we assume to be Gaussian with a ID 
dispersion a. The 2D distribution function of the projected velocity in polar coordinates is 

DF 2 (v 2 )dv 2 d6 = ^ exp(-(6l + v 2 ot - 2v 2 v TOt cos 0)/2)dv 2 d9 , (A3) 

where v 2 = v 2 /a, v TO t = v TO t/cr and 9 is the angle between v 2 and v rot . Upon substitution into 
Eq. [A^ and averaging over all stellar types, one obtains 

— — = 2R s n*a(uoG(v rail x, v m t)) , (A4) 
dr 

where u max = 2uq )S R s / ATa, and 



f»mix I 

G(v msx ,v Iot ) = exp(--D~ ot /2) x y 1 - (x/£ max ) 2 exp(-x /2)I (xv TOt ) dx , (A5) 

where I n is a modified Bessel function of order n. In the simple case where v ro t = 0, 

G(u max ) = ^{) max exp(-^ ax /4)/i (Umax/4) • (A6) 
In the limit of v max — > 00 (AT = 0), 

G(vrot) = \/f exp(-{) r 2 ot /4) [(2 + ?) r 2 ot )/ (i} r 2 ot /4) + 7) r 2 t^i(*rot/4)l , (A7) 
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which is independent of uo :S , so that (uoG(v maiX ,v TOt )) — > («o)G(#rot)- It is straightforward to 
verify that in this limit, the mean transverse velocity approaches its correct asymptotic values, 
since 

lim aG{v rot ) = Jit /2a, 

V rot — >0 V 

lim aG(v Iot ) = v TOt . (A8) 

A fraction of the long duration events are due to stars on trajectories with smaller than average 
impact parameters. Therefore, the median impact parameter of long events is smaller than (uq)/2 
and the median amplification above the threshold is greater than 2. The local average impact 
parameter is 

u(r) = ( / duu DF 2 (v 2 ) d 2 v 2 ) = (u Q W{v m! , x , v TO t)) , (A9) 

J J 

where the weight function W can be written as 

W(v max , v TOt ) = exp(-i) rot /2)- / x(l - (x/v max ) )exp(-x /2)I (xv rot ) dx . (A10) 

2 Jo 

It is straightforward to verify that 

lim W(v mSLX , v rot ) = - . (All) 

fmax-+CiO Z 

The total rate averaged impact parameter is 




(A12) 



and the median amplification above the threshold can be estimated by A = (uq)/u. 
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Figure captions 

Fig. 1. — The GC mass density model. 

Fig. 2. — The KLF of the theoretical continuous star forming model with twice solar metallicity 
(solid line) and solar metallicity (dashed line), compared to the Blum et al. (1996) (3 = 1.875 
power-law model of the composite observed KLF (points). 

Fig. 3. — The differential and integrated microlensing rate as function of the distance behind the 
black hole for the central cluster model, assuming Kq = 16. 5 m and both continuous monitoring and 
AT = 1 yr. Both the total rates and the rates for long duration events with r > AT are shown. 
The discontinuity at 56 pc reflects the transition from unresolved images to the fainter resolved 
images. 

Fig. 4. — The cumulative lensing rate, ri ong (r > AT), and the rate averaged lensing timescale 
^short( T < AT) of the central cluster KLF as function of the sampling interval AT, for the detection 
thresholds Kq = 16 m , 17 m , 18 m and 19 m (from top to bottom). Bold lines show results for faint-star 
lensing and thin dotted lines for bright-star lensing. The thin dashed line in the bottom panels 
is fghort = AT/2. (In the case of resolved lensing, values of f s hort < 0.3 yr are not shown due to 
numerical instabilities in the calculations.) 

Fig. 5. — The mean amplification above the detection threshold, AK, as function of the sampling 
rate AT. Bold lines show results for faint-star lensing and thin dotted lines for bright-star lensing. 
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Figure 5. 



